Dynamical decoupling for a qubit in telegraph-like noises 
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Based on the stochastic theory developed by Kubo and Anderson, we present an exact result of the deco- 
herence function of a qubit in telegraph-like noises under dynamical decoupling control. We prove that for 
telegraph-like noises, the decoherence can be suppressed at most to the third order of the time and the periodic 
Carr-Purcell-Merboom-Gill sequences are the most efficient scheme in protecting the qubit coherence in the 
short-time limit. 
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I. INTRODUCTION 

Dynamical decoupling (DD) is a standard technique to sup- 
press the spin decoherence, with a long tradition in magnetic 
resonance spectroscopy |l|,|2|]. Nowadays, DD is important 
in areas like quantum information processing, since prolong- 
ing the qubit coherence is a fundamental requirement to carry 
out effective operations on quantum states. The basic idea 
of DD is using a sequence of control pulses that frequently 
flip the spins to average out effects of random environmen- 
tal field. DD is originated from Hahn's first spin echo ex- 
periment in 1950 yfl. After Hahn's work, more complex 
pulse sequences were introduced, among which a most fa- 
mous example was the periodic Carr-Purcell-Merboom-Gill 
sequence (CPMG) lillj], which was initially widely used in 
magnetic resonance spectroscopy and, in recent years, was 
introduced to protect the qubit coherence in quantum infor- 
mation processing (3-01 • Concatenated DD (CDD) lH-H] 
is of special interest, since it recursively constructs pulse se- 
quences to eliminate qubit decoherence to an arbitrary order 
of precision in the short-time expansion. The pulse number 
in CDD, however, exponentially increases. The first optimal 
pulse sequence in terms of the pulse number was introduced 
by Uhrig lfl5ll . which was a non-periodic DD. Uhrig's DD 
(UDD) was later shown to be universal in the sense 

that the leading n orders of the time expansion of the decoher- 
ence are eliminated by using n pulses for general finite quan- 
tum systems, i.e., systems with hard high-frequency cutoff in 
noise spectra. Recent experiments achieved remarkable pro- 
gresses in prolonging the coherence via DD schemes IU8I42011 . 

An interesting question is how different DD schemes per- 
form in suppressing spin or qubit decoherence caused by clas- 
sical non-Gaussian noises, such as multi-state telegraph-like 
noises. For Gaussian noises, which are fully characterized 
by their second-order correlation functions, the DD control 
can be readily formulated as the integration of noise spec- 
tra modulated by a filtering function determined by the DD 
sequence J2lll . For general non-Gaussian noises, however, 
such a formalism is not available and one has to rely on cu- 
mulant expansion and/or Gaussian approximation to numeri- 
cally solve the problem. In Ref. 1121 II . numerical calculation 
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shows that CPMG is practically better than CDD and UDD 
for a telegraph noise. Later in Ref. [22], numerical search for 
optimal DD sequences finds solutions close to the CPMG for 
noises with soft cutoffs and to UDD for noises with hard cut- 
offs. These research indicates that CPMG may be optimal in 
combating certain non-Gaussian noises, but more research is 
still needed to reach the conclusive results. 

In this paper, we present an exact result of the decoher- 
ence function under DD control, based on a stochastic the- 
ory (23-25] developed in 1950s mainly by Kubo and Ander- 
son, when they studied the line shape of nuclear magnetic res- 
onance (NMR) spectra. We prove that for a general multi-state 
telegraph-like noise, any DD schemes cannot fully eliminate 
the third order term in the short-time expansion of the deco- 
herence, i.e., the decoherence function is at least of <9(f 3 ), and 
among all possible DD schemes, CPMG is the most optimal 
in suppressing the qubit decoherence in the short-time limit. 
Apart from the theoretical importance, these results are rele- 
vant to the spin decoherence problem in real systems such as 
spins in Si/Si02 interfaces, where the noises from coupling 
with dangling bonds can be approximately represented by a 
discrete multi-state telegraph-like noise 12611 ■ 

The paper is organized as follows. In section|llJ we obtain 
the exact expression for the decoherence function under arbi- 
trary DD control of a qubit in telegraph-like noises. In section 
[Till we expand the decoherence function, solve the optimiza- 
tion problem and prove that CPMG is the global minimum 
solution. 



II. EXACT DECOHERENCE FUNCTION 

We consider a qubit (spin- 1/2) under an external field and a 
random field w(t). The Hamiltonian (in the rotating reference 
frame in which the external field is transformed to be zero) is 

H = S z w(t), (1) 

where S x / y / z is the spin operator along the x/y/z direction. We 
assume w(t) is a multi-state telegraph-like noise, i.e., it jumps 
suddenly and randomly among a set of discrete values {wf), 
with a transition rate matrix T, in the form 

jYj(t) = T jf Y f (f), (2) 

where Yj(t) is the probability for w(f) = wj. Such a noise 
model is widely used in describing various stochastic physical 
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process, such as the spectral diffusion of optical transitions 
due to atom collisions II27I1 . 

The spin coherence is characterized by the transverse po- 
larization x = S x + iS y . The equation of motion for x in the 
form of a so-called Kubo oscillator is 



x(t) = iw(t)x(t). 



(3) 



The ensemble average (x(t)) = J x(t)P[w(t)]D[w(t)] over a 
distribution of the random field /'[w(f)] gives the transverse 
polarization of the spin. Defining Yj(<p, t) as the probability 
density for finding x(f) = e ,(p and w(f) = Wj, we can write the 
ensemble averages (wx) and (x) as 



(wx) 



f)d<p, 



(<p, t)d<p. 



(4a) 
(4b) 



III. OPTIMAL DYNAMICAL DECOUPLING 

To achieve a certain order of DD, we need to find the so- 
lution of the set of (normalized) pulse intervals {a, } to make 
the decoherence function (x(t)) equal to unity up to an error 
of a certain order of t in the short-time expansion. For this 
purpose, we expand the exponential functions in Eq. (O into 
Taylor's series of t. The zeroth order term contains neither W 
norF, and is explicitly 2/>>j(0) = 1 (the sum probability must 
be unity). 

The higher order terms can be classified by the ordering 
of the matrices W and T. A few rules can be established to 
significantly simplify the expansion. First, any term starting 
with F must vanish, because of the probability conservation 
condition 



(9) 



The probability density function satisfies the stochastic Liou- 
ville equation | 



d_ 

di 



Y j {<p,t) = Y J r jj ,Yj,(<p,t) 



(5) 



which contains both the sudden jumps of the random field and 
the precession of the spin polarization. By multiplying e ,ip to 
both sides of Eq. © and integration over the phase angle, we 
obtain 



-y(t) = (T + iW)y(t), 
at 



(6) 



Second, any term ending with F must vanish since 



Ty(0) = 0, 



(10) 



for the random force distribution is stationary. Third, the 
terms of the same order and containing no F sum to zero. This 
is because if F is set to zero (corresponding to the static inho- 
mogeneous broadening condition), the decoherence function 
becomes 



W(a l -a 2 +...+(-l) fi a N+l )t 



i, 



(id 



where yjif) = JC e" p Yj((p,t)dtp, Wjf = SjfWj>. The spin po- 
larization is (x(t)) = Yijyj(t)- 

Under DD control, the spin is subjected to a sequence of flip 
operations between the +z and -z directions. In this paper, for 
the sake of simplicity, we assume the case of ideal ji pulses. 
Equivalently, the flip control can be transformed to the flip of 
the field w(t) in the reference frame rested on the flipped spin. 
Thus, we have the controlled Liouville equation as 



T j(t) = [F + iWf(t)]y(t), 
at 



(7) 



where /(f) is a step-like function jumping between +1 and -1 
every time a control pulse is applied. The exact solution of the 
decoherence function under an -pulse DD control is 



<*( f) > = 2[ e[r+( ~ 



\) N iW]a N +i 



,[r-iW] ai r e ir+iW] 



ait y(0)]., (8) 



where y(0) is the initial probability distribution of the random 
force, and < a„t < 1 is the interval between the nth and 
(n - l)th pulses with 2 a„ = 1. 



under the echo condition a\ — 02 + ... + (-l^a/v+i = 0. 
Using the three rules above, the only non-vanishing term of 
the decoherence up to the third order must have the form of 
Gn{a\, a2, ■ • ■ , a.N)WFWy(Q). To minimize the decoherence 
in the third order, we just need to minimize the coefficient 
G N (a u a 2 ,...,a N ). 

With the conjecture that CPMG could be an optimal solu- 
tion (as in the case of two-pulse control for Gaussian noises 
with hard cutoffs), we write the coefficient Gw(ai, 02, • ■ • , aw) 
as a function of the deviations of the pulse positions from the 
CPMG timing, 



2n-l 
2N 



(12) 



where a n t is the position of the nth control pulse (i.e., a n = 

a„-a„-\). The echo condition is /?i -/?2H \-(-1) n+1 /3n = 0, 

and another constraint is 1 > > a^-i > • ■ ■ > ae\ > 0. 
These two conditions define the physical boundary for an N- 
pulse sequence. Then the coefficient of the third order term 
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+ [ipl [-J3 N + Ifiu-i ~ Wn-i + ... + (-1)%] + 2j3 2 N _i [-ySjv-i + 2fi„- 2 + ■■■ + {-lf~%] 



• + [-A] 2 } 
+ 2yS?(- J 1 )} 



12N 2 



+ h 



N+gN, 



(13) 



where denotes the second-order term, and gN the third- 
order term. The first order term vanishes because of the echo 
condition. Since the second order is always positive (h^ > 0), 
CPMG ifi = 0) must be at least a local minimum solution of 
the decoherence function. 

If CPMG is not the global minimum, we can find a set of 
{fj + 0} so that h N (J3) + g N (J3) < 0. Since h N (J3) > 0, we must 
have gN(fi) < 0. Using a real number A to scale the deviation, 
we get a function of the scaling factor as 

/nW = + h N (A/3) = A 3 g N (/3) + A 2 h N (/3). (14) 

The function /n(A) monotonically decreases with A for A> 1. 
Thus when A is increased from 1, the DD sequence suppresses 
the decoherence better and better. But A cannot be infinitely 
increased under the physical conditions. When A is increased 
to a boundary value Ag, the deviation A/3 will reach the phys- 
ical boundary, at which either two adjacent pulses coincide 
(and become a null operation) or a pulse reaches the boundary 
time at or t. This means a new pulse sequence with fewer 
pulses is obtained. Suppose the new pulse sequence has N' 
pulses (with N' < N), the coefficient can be written as 

1^2 + M*b) = j^- 2 + griff) + h N ,(J3'\ (15) 

where /3' is the deviation from the A^'-pulse CPMG, and g^> 
and htq' are defined as in Eq. (fT3l . Obviously, 

g N ,(J3') + h N ,(J3')<0. (16) 

Then following the same procedure as in the A^-pulse se- 
quence, we can do the scaling A' ft' again to find pulse se- 
quences performing better in suppressing the decoherence by 
increasing A' from 1 to a new boundary value. Then we 
would find a new sequence which performs better with fewer 
pulses than N'. So on and so forth, we will have to con- 
clude that the one-pulse sequence (Hahn echo) is the optimal 
solution among all sequences (including the multiple-pulse 
ones), which can be easily checked to be wrong by comparing 
the performance of the two-pulse CPMG and the Hahn echo. 
Thus, to avoid contradictions, we conclude with the following 
theorem: 

Theorem 1 For an arbitrary multi-state telegraph-like noise, 
CPMG is the global minimum solution of the decoherence 
function among all DD sequences of the same number of 
pulses. 



A corollary which can be directly derived from Eq. (fT3l is 
Corollary 1 The decoherence of a qubit in a multi-state 
telegraph-like noise can not be suppressed by DD beyond the 
third order of the short-time expansion. 

This conclusion is consistent with previous numerical solu- 
tions in Ref. [22], which give pulse sequences close to CPMG 
for boson baths with power-law high-frequency cutoffs. Actu- 
ally, the second order correlation function of the telegraph-like 
noise 

<w(?)w(f')> = J] [We-^WyiO)]. , (17) 

has the form of exponential functions. The Fourier transfor- 
mation to the frequency domain then has a power-law decay 
profile in the high frequency end. Similarly, the higher or- 
der correlation functions (which in general cannot be factor- 
ized into the second order correlation functions as for Gaus- 
sian noises) also have the form of exponential functions. Such 
exponential function form of the correlation functions results 
from the sudden jumps in the telegraph-like noises, which, 
in the physical nature, is induced by instantaneous (Marko- 
vian) collisions in the bath. It is the existence of a very short 
timescale in the system (the collision memory time), or in 
other words, a very large energy scale, that limits the perfor- 
mance of any DD schemes to the third order of precision. 



IV. CONCLUSIONS 

In summary, we have derived an exact expression for 
the decoherence function of a qubit in arbitrary multi-state 
telegraph-like noises, based on the stochastic theory. We 
prove that CPMG is the globally optimal solutions among all 
possible dynamical decoupling sequences of the same number 
of pulses to suppress the decoherence in the short time limit. 
Because of the instantaneous random jumps in the noises, the 
decoherence cannot be eliminated beyond the third order of 
the short time. 
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